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Abstract 

A two-layer statistically equivalent periodic unit cell is offered to predict a macro- 
scopic response of plain weave multilayer carbon-carbon textile composites. Falling- 
short in describing the most typical geometrical imperfections of these material 
systems the original formulation presented in [1] is substantially modified, now al- 
lowing for nesting and mutual shift of individual layers of textile fabric in all three 
directions. Yet, the most valuable asset of the present formulation is seen in the 
possibility of reflecting the influence of negligible meso-scale porosity through a 
system of oblate spheroidal voids introduced in between the two layers of the unit 
cell. Numerical predictions of both the eff'ective thermal conductivities and elas- 
tic stiff'nesses and their comparison with available laboratory data and the results 
derived using the Mori-Tanaka (MT) averaging scheme support credibility of the 
present approach, about as much as the reliability of local mechanical properties 
found from nanoindentation tests performed directly on the analyzed composite 
samples. 

Key words: balanced woven composites, material imperfections, statistically 
equivalent periodic unit cell, image processing. X-ray microtomography, 
nanoindentation, soft computing, numerical homogenization, steady-state heat 
conduction 



^Corresponding author. Tel.: +420-2-2435-4482; fax +420-2-2431-0775 
Email addresses: (Jan Vorel), zeman j @cinl . f sv . cvut . cz (Jan Zeman), 
sejnomOf sv . cvut . cz (Michal Sejnoha), blanka . tomkovaOtul . cz (Blanka Tomkova) 



Preprint submitted to Elsevier 



November 11, 2010 



1. Introduction 

Despite a significant progress in theoretical and computational homogeniza- 
tion methods, material characterization techniques and computational resources, 
the determination of overall response of structural textile composites still remains 
an active research topic in engineering materials science [2]. From a myriad of 
modeling techniques developed in the last decades (see e.g. review papers [3, 4, 
5]), it is generally accepted that detailed discretization techniques, and the Finite 
Element Method (FEM) in particular, remain the most powerful and flexible tools 
available. The major weakness of these methods, however, is the fact that their ac- 
curacy crucially depends on a detailed specification of the complex microstructure 
of a three-dimensional composite, usually based on two-dimensional micrographs 
of material samples, e.g. [6, 7, 8, 5, and reference therein]. Such a step is to a great 
extent complicated by random imperfections resulting from technological opera- 
tions [9, 10], which are difficult to be incorporated to a computational model in a 
well-defined way. If only the overall, or macroscopic, response is the important 
physical variable, it is sufficient to introduce structural imperfections in a cumula- 
tive sense using available averaging schemes such as Voight/Reuss bounds [1 1] or 
the Mori-Tanaka method [12]. When, on the other hand, details of local stress and 
strain fields are required, it is convenient to characterize the mesoscopic material 
heterogeneity by introducing the concept of a Periodic Unit Cell (PUC). 

While application of PUCs in problems of strictly periodic media has a rich 
history, their introduction in the field of random or imperfect microstructures is 
still very much on the frontier, despite the fact that the roots for incorporating ba- 
sic features of random microstructures into the formulation of a PUC were planted 
already in mid 1990s in [13]. Additional extension presented in [14], see also our 
work [15] for an overview, gave then rise to what we now call the concept of 
Statistically Equivalent Periodic Unit Cell (SEPUC). In contrast with traditional 
approaches, where parameters of the unit cell model are directly measured from 
available material samples, the SEPUC approach is based on their statistical char- 
acterization. In particular, the procedure involves three basic steps [15]: 

• To capture the essential features of the heterogeneity pattern, the microstruc- 
ture is characterized using appropriate statistical descriptors. Such data are 
essentially the only input needed for the determination of a unit cell. 

• A geometrical model of a unit cell is formulated and its key parameters are 
postulated. Definition of a suitable unit cell model is a modeling assump- 
tion made by a user, which sets the predictive capacities of SEPUC for an 
analyzed material system. 

• Parameters of the unit cell model are determined by matching the statistics 
of the complex microstructure and an idealized model, respectively. Due to 
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multi-modal character of the objective function, soft-computing global op- 
timization algorithms are usually employed to solve the associated problem. 

It should be emphasized that the introduced concept is strictly based on geo- 
metrical description of random media and as such it is closely related to previous 
works on random media reconstruction, in particular to the Yeong-Torquato algo- 
rithm presented in [16, 17]. Such an approach is fully generic, i.e. independent 
of a physical theory used to model the material response. If needed, additional 
details related to the simulation goals can be incorporated into the procedure with- 
out major difficulties, e.g. [18, 19], but of course at the expense of computational 
complexity and the loss of its generality. 

In the previous work [1], the authors studied the applicability of the SEPUC 
concept for the construction of a single-layer unit cell reflecting selected imper- 
fections typical of textile composites. A detailed numerical studies, based on both 
microstructural criteria and homogenized properties, revealed that while a single- 
ply unit cell can take into account non-uniform layer widths and tow undulation, 
it fails to characterize inter-layer shift and nesting. Here, we propose an extension 
of the original model allowing us to address such imperfections, which have a 
strong influence on the overall response of textile composites [20, 21, 22]. A brief 
summary of the procedure for the determination of the two-ply SEPUC for woven 
composites is given in Section 3. 

Such extensions, however, are hardly sufficient particularly in view of a rel- 
atively high intrinsic porosity of Carbon-Carbon (C/C) composites, which are in 
the center of our current research eff'orts. It has been demonstrated in our previ- 
ous work [23] that unless this subject is properly addressed inadequate results are 
obtained, regardless of how "exact" the geometrical details of the meso- structure 
are represented by the computational model. Unfortunately, the complexity of 
the porous phase seen also in Figure 1 requires some approximations. While 
densely packed transverse cracks aff'ect the homogenized properties of the fiber 
tow through a hierarchical application of the Mori-Tanaka averaging scheme [24], 
large inter-tow vacuoles (crimp voids), attributed to both insufficient impregna- 
tion and thermal treatment, are introduced directly into the originally void-free 
SEPUC in a discrete manner. 

Not only microstructural details but also properties of individual composite 
constituents have a direct impact on the quality of numerical predictions. Infor- 
mation supplied by manufacturers are, however, often insufficient. Moreover, the 
carbon matrix of the composite has properties dependent on particular manufac- 
turing parameters such as the magnitude and durations of the applied tempera- 
ture and pressure. Experimental derivation of some of the parameters is there- 
fore needed. In connection with the elastic properties of the fiber and matrix, the 
nanoindentation tests performed directly on the composite are discussed in Sec- 
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Figure 1 : Representative segment of eight-layer plain weave fabric laminate. 



tion 2 together with the determination of the necessary microstructural parameters 
mentioned already in the previous paragraphs. 

Still, most of the work presented in this paper is computational. In particular, 
a brief summary of the procedure for the determination of the two-ply SEPUC 
for woven composites is given in Section 3. Section 4 is then reserved for the 
validation of the extracted geometrical and material parameters. To that end, the 
heat conduction and classical elasticity homogenization problems are validated 
against available experimental measurements. The concluding remarks and future 
extensions are presented in Section 5. 

2. Experimental program 

As already stated in the introductory part, much of the considered here is pri- 
marily computational. However, no numerical predictions can be certified if not 
supported by proper experimental data [25]. The objective of the experimental 
program in the context of the present study is twofold. First, reliable geometrical 
data for the construction of the unit cell and material parameters of both the car- 
bon fibers and carbon matrix for the prediction of eff'ective properties are needed. 
Since still derived on the basis of various assumptions, these results must be next 
confirmed experimentally to acquire real predictive power. 

Considering the mesoscopic complexity of C/C composites, the supportive 
role of experiments is assumed to have the following four components: 

• Two-dimensional image analysis providing binary bitmaps of the composite 
further exploited in the derivation of two-layer SEPUC 

• X-ray tomography yielding a three-dimensional map of distribution, shape 
and volume fraction of major pores to be introduced into a void-free SEPUC. 

• Nanoindentation tests supplying the local material parameters which either 
depend on the manufacturing process or are not disclosed by the producer. 

For the above purposes a carbon-polymer (C/P) laminated plate was first man- 
ufactured by molding together eight layers of carbon fabric Hexcel G 1169 com- 
posed of carbon multifilament Torayca T 800 HB and impregnated by phenolic 
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resin Umaform LE. A set of twenty specimens having dimensions 25 x 2.5 x 
2.5 mm were then cut out of the laminate and subjected to further treatment (car- 
bonization C at 1000°C, reimpregnation /, recarbonization, second reimpregna- 
tion and final graphitization G at 2200^C (CICICG)) to create the C/C composite, 
see Figure 2 for an illustration and [26] for more details. The reported specimens 
were then fixed into the epoxy resin and subject to curing procedure. In the last 
step, the specimen was subjected to final surface grounding and polishing using 
standard metallographic techniques. 




(a) (b) (c) 

Figure 2: Examples of scanned microstructures; (a) Woven fabric, (b) carbonized composite, 
(c) graphitized composite 



2.7. Two-dimensional image analysis 

The actual image analysis device used for structural image acquisition and 
analysis consists of NIKON ECLIPSE E 600 microscope, Marzhauser motorized 
scanning stage, digital monochrome camera VDC 1300C and image analysis soft- 
ware LUCIA Note that high reflectance of the woven fabric allows relatively 
good visual resolution of individual parts of a composite structure as demonstrated 
in Figure 2. Using the method of gradual abrasion, transverse sections of the 
composite laminate allowed us to generate a database of micro-images intended 
for further processing. Unfortunately, a low color contrast of the reinforcement 
(carbon fabric) and matrix is a major impediment to an automatic separation of 
individual objects. Therefore, a manual preprocessing of images by marking the 
borders of selected objects becomes often necessary, cf. Figure 3. Further image 
analysis and object measurement was, however, fully automatic providing an av- 
erage thickness of carbon tows, shape and dimensions of fiber tow cross-section, 
size and shape of major voids, distribution of transverse and delamination cracks, 
etc. 

Although, as demonstrated in Section 4, statistical consideration of the ac- 
quired results proved useful, the essential input for the preparation of physical 
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(a) 



(b) 



Figure 3: Image processing steps; (a) Manual marking of object borders, (b) automatic marking 
of objects 

models for computational analysis (SEPUC) is given by binary images of the com- 
posite. An illustrative example of the result of two-dimensional image processing 
is available in Figure 4 implicating that material porosity is neglected when treat- 
ing only the geometrical imperfections of the fabric reinforcements. At this point, 
a direct comparison of the resulting mesoscopic predictions with experimental 
measurements would therefore be meaningless. 




(b) 

Figure 4: Examples of image analysis; (a) Woven fabric, (b) multi-layered C/C composite 

2.2. Three-dimensional X-ray microtomography 

Porosity of C/C composites plays an important role in the derivation of ef- 
fective material properties [27]. Most common approach to characterizing this 
property employs sectioning, recall Section 2.1. The influence of shape of pores, 
estimated from 2D images of real C/C composites [28], on the mechanical re- 
sponse has been addressed e.g. in [29, 27]. The X-Ray microtomography becomes 
a valuable tool rendering three-dimensional phase information [30, 31, 32]. In the 
present study, high resolution computer tomography images provided by the Inter- 
facultary Laboratory for Micro- and Nanomechanics of Biological and Biomimet- 
ical Materials of the Institute of Lightweight Design and Structural Biomechanics 
were used to obtain the shape, size, location and volume fraction of inter-layer 
(crimp) voids. A particular example of the distribution of major porosity in C/C 
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multi-layered plane- weave composite is presented in Figure 5. While the basic 
characteristics of the porosity can be directly extracted from these images, a direct 
introduction of pores in their full complexity as seen in Figure 5(a) is impossible. 
Instead, a set of oblate spheroids approximating the shape and volume of actual 
pores is accommodated in between the two layers of the SEPUC fitting their true 
location as close as possible. The results of two- and three-dimensional image 
analysis appear in Table 1, quantifying the structural parameters of a sigle-layer 
textile composite, mutual shift of individual layers and the volume and shape of 
macroscopic pores. 




(a) (b) 

Figure 5: X-ray microtomography; (a) Interior distribution and shape of large vacuoles, (b) three- 
dimensional view of the porous composite structure 



Parameter 




Average 


Standard deviation 


Tow period 


[/im] 


4500 


300 


Tow height 


[/um] 


150 


20 


Inter-tow gap 


[jum] 


400 


105 


Layer height 


[jum] 


300 


50 


Horizontal shift 


[jum] 





675 


Vertical shift 


[/um] 





110 


Porosity 


[%] 


8 


3.5 


Pore aspect ratio 


[-] 


0.4 


0.2 



Table 1 : Parameters of the periodic unit cell [28] 



2.3. Phase elastic moduli from nanoindentation 

Prediction of complex macroscopic response of highly heterogeneous mate- 
rials from local phase constitutive theories is a an important aspect of microme- 
chanical modeling. The reliability of these predictions, however, is considerably 
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influenced by available information on material data of individual constituents. 
Even though supplied by the producer, these information are often insufficient for 
three-dimensional analysis. It is also known that material properties of the matrix 
much depend on the fabrication of composite and may considerably deviate from 
those found experimentally for large unconstrained material samples [33]. 

Carbon matrix developed in the composite through a repeated process of im- 
pregnation, curing and carbonization of the phenolic resin is a solid example. This 
resin belongs to the non-graphitizing resins so that the final carbon matrix essen- 
tially complies, at least in terms of its structure, with the original cross-linked 
polymeric precursor. Therefore, the resulting material symmetry is more or less 
isotropic with material parameters corresponding to those of glassy carbon. Nev- 
ertheless, when constrained the assumed matrix isotropy may evolve into the one 
of the fibers particularly in their vicinity. Although the PAN (polyacrylonitril) 
based carbon fibers are known to have a relatively low orderliness of graphen 
planes on nano-scale, they still posses a transverse isotropy with the value of lon- 
gitudinal tensile modulus (usually available) considerably exceeding the one in 
the transverse direction (often lacking). Additional experiments, preferably per- 
formed directly on the composite, are therefore often needed to either validate the 
available local data or to derive the missing ones. 

At present, nanoindentation is the only experimental technique that can be 
used for direct measurement of mechanical properties at material micro-level. A 
successful application of nanoindentation to C/C composites has been reported 
in [34, 35, 36]. In the present study, our attention was limited to the evaluation 
of the matrix elastic modulus and the transverse elastic modulus of the fiber. The 
remaining data were estimated from those available in the literature for similar 
material systems. 




(a) (b) 

Figure 6: Nanoindentation - location of indents; (a) Transverse direction, (b) longitudinal direc- 
tion (compression) 



Three locations, distinctly separated in optic microscope, were tested - ma- 
trix, parallel fibers Figure 6(a), perpendicular fibers Figure 6(b). The matrix was 
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therefore assumed isotropic and possible anisotropy, which may arise inside the 
fiber tow, was not considered. As seen in Figure 6 several measurements were 
recorded for each of the three locations. The measurements were performed us- 
ing CSM Nanohardness tester equipped with a Berkovich tip allowing for 0.1-500 
mN loading range. To ensure elastic response relatively low indentation forces up 
to lOmN were considered. The elastic moduli were extracted from an unloading 
part of the indentation curve using the well known Oliver-Pharr procedure [37]. 
In particular, the indentation elastic modulus is then provided by 

where A is the projected contact area at the peak load and S is the contact stiff'ness 
evaluated as the initial slope of the unloading curve. The following equation is 
then adopted to account for a finite elastic stiff'ness of the indenter 

1 1 - y2 1 - 
— = ^ + ^ , 

tested material and Ei and Vi are parameters of the indenter (for diamond: Ei=ll4l 
GPa and =0.07). In this study, the matrix Poisson ratio was assumed equal to 0.2 
while the fiber Poisson ration was set equal to 0.4. The complete set of parameters, 
both measured averages labeled by * and those adopted from the literature, is 
available in Table 2. Note that the matrix modulus agrees relatively well with the 
one found for the glassy carbon in [36]. 





Young modulus 


Shear modulus 


Poisson ratio 


Material 


[GPa] 


[GPa] 


[-] 


fiber longitudinal 


294 


11.8 


0.24 


transverse 


12.8* 


4.6 


0.4 


matrix 


23.6* 


9.8 


0.2 



Table 2: Material parameters of individual phases 



The final note is concerned with heat treatment of the C/C composite dur- 
ing fabrication. It has been observed experimentally [38] that even for T800 
based composites the tensile Young's modulus increases for graphitized speci- 
mens (CICICG) when compared to only carbonized ones (CICIC). It is suggested 
that this phenomenon may be caused by further stiff'ening of carbon fibers. This, 
however, is difficult to address in the present study as the tensile modulus cannot 
be measured via nanoindentation. On the contrary, neither the matrix properties 
nor the porosity profile is expected to change considerably with graphitization. 
Therefore, the experimental data reported in [38] for the CICIC system rather 
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than those for CICICG system will be fostered for comparison with numerical 
predictions. 

2.4. Laboratory evaluation of effective properties 

Corroboration of a mechanics model by experimental data is still thought vital 
for the model to be accepted, inasmuch as there is simply nothing better, even 
though an experiment often comprises laboratory measurements and a theory for 
calculating not directly measurable parameters. 

With regard to thermophysical parameters the pulse transient method [39] 
combined with a heat loss model for the calculation of temperature response, 
when lower currents are used for pulse generation, is adopted. The searched 
thermal conductivities are found subsequently using a general relation between 
thermal conductivity, specific heat and thermal difflisivity, the latter two given in 
terms of the maximum attained temperature, the elapsed time to reach this temper- 
ature, total amount of generated heat and a set of correction factors characterizing 
deviations of calculated thermal diff'usivity and specific heat from those derived 
for an ideal system (ideal heat source placed in bundleless isotropic body). It is 
not the objective of this section to provide all details of this method, as these can 
be found, e.g. in [39, 40], but rather to suggest the complexity of deriving the 
eff'ective thermal conductivities from an experiment, which in turn are to be used 
to question the quality of numerical predictions presented later in Section 4. In 
particular, the results presented in [40] (C/C laminate (EXP) row) are used to val- 
idate the numerical model. Averages of the measurements from six samples of an 
eight layer C/C laminate are available in Table 3^. 



material Thermal conductivities [Wm ^] 





Xn 


X22 


/f33 


air 


0.02 


0.02 


0.02 


fiber 


35 


0.35 


0.35 


matrix 


6.3 


6.3 


6.3 


porous tow (MT) 


24.12 


1.05 


1.42 


C/C laminate (EXP) 


10 (warp) 


10 (weft) 


1.6 



Table 3: Phase [26], unidirectional (UD) C/C composite (porous fiber tow) and laminate eff'ective 
thermal conductivities [40] 

Similar difficulties arise when deriving the dynamic tensile and shear moduli 
from a resonant frequency method. This technique was employed in [38] to derive 



^See Figure 7 for the definition of local x (fiber tow) and global X (textile ply) coordinate 
system. Note that the local xi axis is aligned with the fiber direction and the global in-plane Xi 
and X2 axes represent warp and weft directions, respectively. 
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the homogenized tensile and shear moduh of a four layer T800 fiber fabric based 
plain weave C/C composite. As part of their study a unidirectional (UD) car- 
bonized CICIC composite was examined. The resulting moduli stored in Table 3 
(porous tow (EXP) column) were utilized here to validate the Mori-Tanaka esti- 
mates of the homogenized properties of the porous fiber tow (porous tow (MT) 
column). Note also the corresponding MT predictions for the heat conduction 
problem listed in Table 3. 

Contemplation of the MT micromechanics model on the level of fiber tow has 
two reasons. First, it provides the homogenized fiber tow conductivity and stiff'- 
ness matrices needed in the mesoscopic study in Section 4 relatively easily. Sec- 
ond, the acquired predictions are sufficiently "reliable" when no microstructure 
information other than volume fractions of individual constituents are at hand. In 
the present study a relatively high volume fractions of fibers equal to 63.5% and 
open porosity by water penetration of 6.6% were taken from [38]. Comparing the 
predicted and measured longitudinal moduli in first two columns in Table 4 both 
verifies the local phase properties given in Table 2 and justifies the use of the MT 
method. It goes beyond the present scope, however, to give full details on this 
method. A few remarks are presented in Section 4.3, further details are available 
in [41, 12, 24, to cite a few]. 



parameter 


porous tow (MT) 


porous tow (EXP) 


C/C laminate (EXP) 


En 


193.8 


«200 




Gi2 


10.3 


« 11.5 






8.0 






E33 


14.3 






Gi3 


7.4 






G23 


4.0 






V12 


0.23 






Vl3 


0.23 






V23 


0.38 







Table 4: Effective elastic properties of UD C/C composite (porous fiber tow) and C/C laminate 
conductivities [38]. Young's moduli are given in [GPa] 



3. Statistically equivalent period unit cell 

The concept of Statistically Equivalent Periodic Unit Cell for random or im- 
perfect microstructures is now well established. Individual steps, enabling the 
substitution of real microstructures by their simplified artificial representatives - 
the SEPUCs - are described, e.g. in [14, 1, 15] and additional references given 



11 



below. Herein, these steps are briefly reviewed concentrating on the specifics of 
multi-layer woven composites. 

3.1. Geometrical mesostructural model 

The basic building block of the adopted SEPUC is provided by a single-ply 
model of plain weave composite geometry proposed by Kuhn and Charalambides 
in [42]. The model consists of two orthogonal warp and weft tows embedded in 
the matrix phase and it is parametrized by four basic quantities, directly measur- 
able by two-dimensional image analysis (recall Table 1): the half-period of tow 
undulation a, the maximal tow thickness b, the width of the intra-tow gap g and 
the overall height of the ply h, see Figure 7 (a). The three-dimensional woven 
composite SEPUC, shown in Figure 7 (b), is formed by two identical one-layer 
blocks, relatively shifted by distances Ai, A2 and A3 in the direction of the cor- 
responding coordinate axes. Finally, cutting a SEPUC by the plane X2 = a or 
Xi = a yields the warp or weft two-dimensional sections, used as the basis for the 
determination of the unit cell parameters. 




(c) 

Figure 7: Geometrical model of SEPUC; (a) Two-dimensional cut of a one-layer model, (b) two- 
layer model including periodic extension of upper layer, (c) two-dimensional cut 



12 



3.2. Quantification of random microstructure 

Assuming a statistically homogeneous and ergodic binary material system, 
two basic statistical functions are available to capture essential characteristics of 
the analyzed tow-matrix material system. The first descriptor is a two-point proba- 
bility function S {X) [43], which quantifies the probability of two points, separated 
by a vector Z, being both found in the domain occupied by the warp and weft tows. 
The alternative statistics, proposed by Lu and Torquato [44] to capture long-range 
eff'ects, is the linear path function L{X) giving the probability that a randomly 
placed segment X is fully contained in the tow region. Both descriptors can 
be easily computed for digitized microstructures; in particular, the Fast Fourier 
transform library FFTW [45] is used to evaluate the S function and the sampling 
template consisting of Nd concentric rays discretized by A^^ pixels (cf. Figure 8) is 
employed to determine the linear path function. The periodic boundary conditions 
were adopted for both descriptors to eliminate edge eff'ects [46]. 



Figure 8: Example of sample template for lineal path function generated by A^^ = 16 segments, 
each of which is discretized using = 11 pixels. 



3.3. Calibration ofSEPUC parameters 

In overall, the adopted model of the unit cell involves seven independent pa- 
rameters 

y = [a,Z?,g,/i,Ai,A2,A3], (1) 

to be determined from available microstructural data. For the sake of generality, 
we assume that the microstructure configuration is characterized by microstruc- 
tural function associated with (at most) warp and weft directions; i.e. functions 
5'warp and Lwarp for the warp cross-section and descriptors 5 weft and L^eft for the 
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weft cross-section, recall Figure 7(b). ^ In particular, see also [15, 16], the follow- 
ing quantities are introduced to measure the similarity between a SEPUC and the 
original microstructure: 



Fjiy) = 



1 



E E Z (Sp(y^iJ)-Sp(iJ)) , (2) 



•maxjmax 



p€{warp,Weft) (=-imax i=-7max 



Fliy) = 



NdNt 



J]J]{L,(y,i,j)-L,{i,j)) , (3) 



pe{warp,weft} i=l j=l 



where, e.g. 5'warp(y, ^ j) denotes the two-point probability function determined 
for the warp cross-section of a SEPUC described by parameters y and the value 
of argument X = [/, j], LweftCy, ^ j) stores the value of the weft-section lineal path 
function for the j-th pixel of the i-th segment, S. and L. denote the statistics 
related to original media and the dimensions /max and jmax are determined as half 
of the minimum height and width of the bitmaps representing a SEPUC and the 
reference image. 

The two-dimensional data, obtained via image analysis in Section 2.1, can be 
complemented by independent experimental measurements of three-dimensional 
volume fractions of the tow phase, see [26] for further details. Such information 
is accounted for by an additional discrepancy measure 



where (p(y) denotes the SEPUC three-dimensional volume fraction and (p is the tar- 
get value. The former quantity is determined from the analytical representation of 
the SEPUC geometry [42] using an adaptive Simpson quadrature [47, Chapter 4] 
with the relative accuracy set to 10"^. 

Moreover, the multiple descriptors can be arbitrary combined in the form of a 
weighted sum. For example, if all available information is employed, the objective 
function attains the form 



with a. denoting scale factors used to normalize the influence of each descriptor, 
determined from twenty randomly generated SEPUCs in the current study. 

The final term is introduced into the objective function to eliminate the in- 
tersection of the upper-layer and lower-layer tows. To that end, we compute the 



^When the same statistics is assumed for both warp and weft directions we set Ai = A2. 



(4) 



(5) 
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overlap 6 as the minimum signed distance between the upper and lower tow sur- 
faces and introduce the constraint 5 > via a polynomial exterior penalty: 

/D(j) = (l + ^j^/^D(j), (6) 

where 5_ denotes the negative part of 5, D refers to a particular combination of the 
descriptors and the value of exponent is set to /? = 3. Note that this approach was 
inspired by a recent work of Collins et al. [48] related to high-density polydisperse 
particulate composites. 

Now, the optimal values of the SEPUC parameters can be determined as the 
solution to a box-constrained global optimization problem 

y G Argmin/£>(z), (7) 
l<z<u 

where the lower and upper bounds / and u are directly based on the image anal- 
ysis data acquired in Section 2.1. A closer inspection reveals that objective func- 
tions (6) are multi-modal and discontinuous due to the effect of limited bitmap 
resolution, cf. [49]. Based on our previous experience with evolutionary optimiza- 
tion, a stochastic optimization algorithm RASA [50, 51], based on a combination 
of a real- valued genetic algorithm and the Simulated Annealing method, is used 
to solve the optimization problem (7). 

3.4. Verification of optimization procedure 

Before applying the proposed methodology to the multi-layered C/C system, 
we first investigate the robustness of the stochastic optimization algorithm and the 
accuracy of the identified parameters. To that end, an artificial two-layer compos- 
ite, with structural parameters corresponding to the average values in Table 1 (i.e. 
a = 2, 250 yum, Z? = 150 yum, g - 400 yum and h - 300 yum) and the layer shifts 
Al = A2 = 1, 125 yum and A3 = -80 yum, is selected as the target material. The unit 
cell cross-section is digitized into a 1, 024 x 140 pixel bitmap, which corresponds 
to horizontal and vertical resolutions of ^ 4.4 yum per pixel. A template with pa- 
rameters A^j = 16 and A^^ = 70 pixels is used to sample the values of the lineal 
path function. As the focus of this section is on the influence of the statistical de- 
scriptor on the accuracy of SEPUC parameters determined from planar statistics, 
only two-point probability and lineal path functions are employed in the current 
example. 

Due to the stochastic nature of the optimization algorithm, we present the 
statistics of the results obtained from twenty independent executions. For each 
component of the vector j, the lower and upper bounds were set to 50% and 
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200% of the target value, respectively. The algorithm was terminated once a solu- 
tion with the value fo smaller then 10"^ was found. The setting of all parameters 
of the RASA algorithm are specified in detail in [49, page 65]. The results appear 
in Tables 5 and 6, storing average values and standard deviations of the number 
of the objective functions evaluations and of the geometrical parameters, respec- 
tively. Note that the standard deviations arise due to stopping criteria and due to 
the limited bitmap resolution. 



Objective 


Success rate 


Number of function calls 






Average Standard deviation 


fs 


20/20 


14,586 13,115 


ft 


20/20 


3,698 1,231 


fs+L 


20/20 


9,446 1,614 



Table 5: Verification example: Success rate and statistics of the number of objective function 
evaluations 



Objective 


a [//m] 


b [yum] 


g [yum] 


h [/um] 


Ai = ^2 [/^m] 


A3 [//m] 


fs 


2,245.9 


149.7 


394.5 


299.8 


1, 125.9 


-79.1 




(6.0) 


(0.2) 


(7.9) 


(0.6) 


(9.2) 


(1.8) 


ft 


2, 249.6 


151.8 


368.5 


311.0 


1,184.2 


-107.0 




(44.4) 


(1.4) 


(37.5) 


(14.5) 


(152.5) 


(34.4) 


fs+L 


2, 249.9 


150.0 


399.9 


302.7 


1,127.2 


-82.0 




(6.8) 


(0.1) 


(4.1) 


(2.0) 


(4.4) 


(3.3) 



Table 6: Verification example: Accuracy of identified parameters; numbers in parentheses corre- 
spond to standard deviations 



For the two-point probability-based objective function fs, we observe the pa- 
rameters of the artificial target unit cell are back-identified with accuracy superior 
to the bitmap resolution 4.4 //m per pixel. From all the SEPUC parameters, the 
inter-tow gap g and the horizontal shift Ai = A2 display the highest scatter; its 
value is, however, comparable again to the pixel size. The major difficulty with 
the fs objective function is a large scatter in the required number of function eval- 
uations, cf. Table 5, which is the consequence of its highly multi-modal character. 
When the lineal path function is employed as the only objective, the optimization 
problem becomes easier to be solved; the average number of iterations consid- 
erably decreases to approximately 25% and the coefficient of variation reduces 
by about 60%, but it comes at the expense of substantially higher scatter in the 
optimal SEPUC parameters. Similarly to, e.g. [17, 52, 53], we observe the com- 
bination of both statistics results in pixel accuracy and scatter in all the parame- 
ters for a moderate number of objective function evaluations and therefore both 
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functions will be employed for to determine SEPUC parameters for disordered 
micrographs. It is also worth noting that, similarly to [1], the global optimization 
algorithm succeeded in locating the global optimum for every execution of the 
algorithm. This demonstrates the robustness of the RASA optimizer and further 
justifies the application of heuristic algorithms to the SEPUC determination. 

3.5. SEPUC for multi-layered C/C composite 

Having verified the optimization algorithm for the SEPUC calibration, we now 
proceed with the analysis of the eight-layer C/C composite represented by bitmap 
appearing in Figure 4(c). First, to keep the optimization process manageable, the 
original 2, 261 x 861 bitmap was down-sampled to a 1, 024 x 390 image, resulting 
again in the pixel size of 4.4 yum. In accordance with conclusions of the previous 
section, all available information was employed to determine SEPUC parameters, 
constrained to the range ± 3cr/, with y- and cr^ denoting the mean and the standard 
deviation of the /-th structural parameter taken from Table 1 . Parameters of the 
algorithm were set to the identical values as previously and the termination cri- 
terion was set to 50,000 of fs+L+(p function evaluations. Again, the optimization 
was executed independently twenty times to obtain reliable results. 

The resulting two-point and lineal-path functions corresponding to the SEPUC 
and the original microstructure appear in Figure 9. In the Xi-axis direction, the 
original statistics is well-reproduced, particularly at the X3 = //m plane where 
the extreme values and the shape of the descriptors are in almost perfect agree- 
ment. The SEPUC also partially captures the local peaks appearing in the multi- 
layer system for ^ +200 yum, cf. Figure 9 (a). In the perpendicular direction, 
we observe that the SEPUC two-point probability function is influenced by pe- 
riodic boundary conditions, leading to more oscillatory behavior and to a slight 
shift of the minimum values from 41% to 35%. The match between the lineal 
path functions in terms of extreme values and shape of the function is even closer, 
since this descriptor is non-periodic even for periodic microstructures. Finally 
note that the tree-dimensional volume fraction (p of the SEPUC is 51%, which co- 
incides exactly with the reference value taken from [26]. Therefore, we conjecture 
that the idealized SEPUC captures the dominant geometrical features of the orig- 
inal system; the diff'erences visible from Figure 9 arise mainly due to the periodic 
boundary conditions and idealized shape of SEPUC (see also Figure 10(b)). 

In Table 7, we present the parameters of the SEPUC together with the standard 
deviations estimated from independent optimization runs. The scatter in the iden- 
tified parameters is comparable to the values reported in the verification example 
and supports the claim that the SEPUC corresponds to the global minimum of 
the optimization problem. The values of the parameters demonstrate that SEPUC 
captures a moderate horizontal shift of individual layers and their mutual overlap. 
These conclusions are further supported by a three-dimensional representation 
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(a) (b) 




(c) 

Figure 9: Statistical descriptors for multi-layered C/C composite and SEPUC; (a) two-point prob- 
ability function and (b) lineal path function determined for multi-layered composite, (c) two-point 
probability function and (d) lineal path function corresponding to SEPUC 



in Figure 10(b), showing that SEPUC reproduces the matrix rich regions together 
with the strong nesting of individual layers, of course within the constraints of the 
selected geometrical model and the tow impenetrability condition. 

3.6. Computational model 

The end goal of Section 3 is the formulation of a computational model in- 
tended for the finite element based homogenization. This implies the use of con- 
forming finite element meshes enabling the implementation of periodic boundary 
conditions mentioned later in Section 4.1. This might seem daunting in that it 
requires not only incorporation of an arbitrary shift of the two layers of fabric 
reinforcement, but also an independent introduction of voids discussed in Sec- 
tion 2.2. An illustrative example of the geometry of such a model is presented 
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a [/um] 


b [/um] 


g [/^m] 


h [yum] 


Ai = ^2 [/im] 


A3 [/mi] 


2,181 


118 


394 


251 


288 


-47 


(10.0) 


(0.3) 


(0.1) 


(1.3) 


(10.5) 


(2.0) 



Table 7: Optimal parameters of the two-layer periodic unit cell free of pores; numbers in paren- 
theses correspond to standard deviations 



in Figure 10. 




(a) (b) 

Figure 10: Computational model; (a) Two-dimensional cut of a two-layer model with voids (b) 3D 
view of the geometry of a two-layer UC model with voids 



In the present study this step is accomplished by combining the principles 
of CAD modeling [7] with the volumetric modeling capacities of the ANSYS® 
package"^. 

In order to ensure the symmetry of the resulting FEM mesh, a primitive block 
of the tow is modeled first, see Figure 1 1(a). Subsequently, using mirroring, copy- 
ing and merging operations, the whole volume of one reinforcement layer is gen- 
erated. Figure 11(b). The second layer is derived analogously and then placed 
according to parameters Ai,A2 and A3 as shown in Figure 11(c). The porous 
phase is introduced next being represented by four identical oblate spheroids, the 
volume of which is derived from X-ray microimages shown in Figure 5. These 
are then periodically extended over the entire model. Figure 1 1(d). Their location 
is assumed to mimic the distribution of large vacuoles that typically appear, as 
also seen in Figure 1, in the location of tow crossings. However, this is difficult 
to achieve in general, and this is also why the porous phase was excluded from 
the minimization problem presented in Section 3. Finally, the volume correspond- 
ing to the matrix is generated using the subtraction of the body of reinforcements 
and the 2a x 2a x (2h + A3) parallelepiped defining the SEPUC. The resulting 
geometrical model appears in Figure 10(b). 



www. ansys . com 
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Now, the mapped meshing technique [6, 54] can be employed to ensure pe- 
riodicity of the resulting finite element mesh. First, half of the external surfaces 
of the SEPUC are discretized and the mesh is then copied to the homologous sur- 
faces. Next, the tetrahedral elements corresponding to tows, voids and matrix are 
generated based on the data created in the previous steps. The corresponding finite 
element mesh is shown for illustration in Figs. 1 l(e)-(f). 




(e) (f) 

Figure 1 1 : Finite element mesh generation; (a) primitive volume, (b) one layer of reinforcements, 
(c) two layers of reinforcements, (d) two layers of reinforcements with voids, (e) FEM mesh of 
tows and voids, (f) FEM mesh of SEPUC 



4. Numerical evaluation of effective properties 

Numerical evaluation of eff'ective elastic moduli and thermal conductivities, 
the most classical subject in micromechanics, is descried in this section. This se- 
lection of mechanical and heat conduction problem is promoted not only by avail- 
able experimental measurements but also by their formal similarity, considerably 
simplifying the theoretical treatment as seen hereinafter. 
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4.1. Theoretical formulation of homogenization 

First-order homogenization approaches are now well established and described 
in many papers [55, 56, to cite a few] to provide estimates of effective properties 
of material systems with periodic microstructures. Bearing in mind the analogy 
between basic quantities related to heat conduction problems (i.e. the local mi- 
croscopic h(x) and uniform macroscopic H temperature gradients and fluxes q(x) 
and Q as their conjugate measures) and the corresponding quantities applied to 
mechanical problems (local s(x) and uniform E strains and the conjugate stress 
measures cr(x) and E), we consider a heterogeneous periodic unit cell Y and vari- 
ations of local temperature 6(X) and displacement u(X) fields written in terms of 
the uniform macroscopic quantities H and E as 

e(X) = H 6\Xl u(X) = E X^ w*(Z), (8) 

where 6* and u* are F-periodic temperature and displacements fluctuations, re- 
spectively and •(X) is introduced to represent a given quantity in the global co- 
ordinate system Z, recall Figure 7. Next, denoting xi^) the local conductivity 
matrix and similarly L(x) the local stiff'ness matrix, the local microscopic consti- 
tutive equations in the local coordinate system x become 

= -x(x)h(x), o-(x) = L(x)s(x). (9) 

To complete the set of equations needed in the derivation of eff'ective properties, 
we recall the Hill lemma for mechanical problem and the Fourier inequality [57] 
for an equivalent representation of steady state heat conduction problem together 
with Eqs. (9) and write the global-local variational principles, see e.g. [23, 58, for 
further details] in the forms 

{6h(xyx(xMx)) = 0, {6s(xyUx)s(x)) = 0, (10) 

where {a(x)) represents the volume average of a given quantity, i.e. {a(x)) = 
a(x) dQ. In the framework of finite element approximation, the discrete 
forms of local gradients derived from Eqs. (8) read 

h(X) = H + B^(Z)6/* , s(X) = E + B"(Z)w* , (1 1) 

where B* stores the derivatives of the element shape functions w.r.t. X and 6*^ and 
w* are the vectors of the fluctuation part of nodal temperatures and displacements, 
respectively. Substituting Eqs. (11) into Eqs. (10) gives 

6e:;'{BWx(X)B\x))e: = -j^f (b^(Z)V^)}h, (12) 

(b"(Z)^L(Z)B"(Z)}w* = (b"(Z)^L(Z)}£, (13) 
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to be solved for nodal temperatures 0^ and nodal displacements w*. Combining 
Eqs. (11) and Eqs. (9) now allows us to write the volume averages of local heat 
fluxes and local stresses as 

Q = (TWq(x)) = -7^ r T\xyx(x)T\X)h(X)dQ, (14) 

E = (t^(Z)V(x)) = t^. f T^(X)'Ux)T^(X)e(X)dQ, (15) 

also showing the relationship between material matrices in the local and global 
coordinate systems in terms of transformation matrices T^, T^, see e.g. [59, 49, 
23, 24], 

X(X) = TWx(x)T\XI L(Z) = T^(Z)TL(x)T^(Z). (16) 

The results of Eqs. (14) and (15) renders the macroscopic constitutive laws in the 
form 

Q = -x^'H, E = L^E, (17) 

where x^ ^i^d L"^ are the searched homogenized effective thermal conductivity 
and elastic stiffness matrices, respectively. In particular, for a three-dimensional 
SEPUC the components of the 3 x 3 conductivity matrix follow directly from 
the solution of three successive steady state heat conduction problems. To that 
end, the periodic unit cell is loaded, in turn, by each of the three components of 
H, while the other two vanish. The volume flux averages, Eq. (14), normalized 
with respect to H then furnish individual columns ofx^- The components of the 
6x6 elastic stiffness matrix L"^ are found analogously from the solution of six 
independent elasticity problems together with Eq. (15). 

4.2. Numerical simulations 

Section 3.6 introduced the computational model of a two-layer SEPUC on 
the basis of geometrical parameters derived in Section 3.5, see also Table 7. 
The mesoscopic porosity (large intertow vacuoles seen in Figs. 1 and 5) equal 
to 5.5% was subdivided, as suggested in Section 3.6, into four oblate spheroids 
(^1 = ^2 > ^3) given the ratio of principal semi-axes as high as possible to fit in 
between the two layers. Note that this value corresponds to the volume of pores 
derived for the four-ply laminate in Figure 1 to make these results comparable 
with those presented in the next section, where the same laminate is analyzed us- 
ing the Mori-Tanaka method. It is fair to mention that this value is slightly lower 
than the one provided by X-ray microtomography, which, on the other hand, is 
surely enhanced by an "artificial" open porosity in the vicinity of specimen edges 
caused by mechanical handling of the specimen before taking the measurements. 
Apart from that, the X-ray microimages still show a shape of pores similar to an 
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oblate spheroid and their ahgnment with layers of carbon tows as assumed in the 
present study. 

Having the computational model the actual derivation of effective properties 
follows the steps outlined in the last paragraph of the previous section. The re- 
sults pertinent to both heat conduction and elasticity problems are summarized 
in Table 8. 



Thermal conductivities 


Elastic moduli 


[Wm-iK-i] 


[GPa] 


/l^ warp/ weft A^trans 


£^11 G\2 £'33 G\3 


9.0 1.9 


60.3 8.1 11.9 5.3 



Table 8: The homogenized thermal conductivities and elastic moduli derived from the application 
of the two-layer SEPUC plotted in Figure 1 1(f) and FEM. 



4.3. Mori-Tanaka approximations 

The Mori-Tanaka method is known to provide quick and reliable estimates 
of the effective properties of heterogeneous materials albeit having only limited 
information about their microstructure. Although extensively used for several 
decades, it was only recently when this method was shown to provide reason- 
able results even for imperfect highly porous textile C/C composites. Inasmuch 
this method is used in the present study to supply, apart from laboratory measure- 
ments, additional results to compare with FEM simulations, we take the liberty 
of giving a brief summary of this method on the subject of meso-scale problem 
in Fig. 7 and refer the interested reader to [60, 12, 24]. We begin by writing the 
homogenized conductivity matrix of C/C composites, first in the absence of the 
porous phase, in terms of the partial temperature concentration factors (localiza- 

r^,,p(0,(A,0) and weft T^^,(f, 



tion tensors) for the warp Tl^^^JO, if/, 0) and weft . (|, if/, 0) directions as 



Xm X\ ~^ 2 {[((/^warp'^warp)) X\ ((^warp))] ~'~ 

+ [({;twefXeft»-Ari({TL«))l)x 

X (cil+f [({TLp)) + {(ne«))]}"\ (18) 

where is the 3 x3 conductivity matrix of the matrix phase and the double brack- 
ets (( )) denote averaging over all possible orientations depending on the adopted 
orientation distribution function if/, ^) with 0, if/ and ^ being the Euler angles. 
This ushers in the issue of random distribution of orientation angle if/ (see also 
Fig. 7(a)) along the fiber tow path addressed already in [12, 24]. Providing the 



23 



X 



histograms of inclination angle if/ are available, consult e.g. [12] for that matter, 
the orientation average of, say Twarp(0, ij/u 0), then reads 

m 

((Twarp)) = Yj ATwarp(0, 0), (19) 

where m denotes the number of sampling values. The discrete angles ij/i and prob- 
abilities Pi follow directly from the image analysis data [61]. 

Formal similarity between heat conduction and elasticity problems then pro- 
vides the homogenized stiffness matrix in the form, recall Eq. (18), 

= y {[((LwarpT^arp)) " Li ((X^arp))] + 

+ [((LweftT^eft)) - Ll ((T^eft))]) X 

[c^^^'ii^^^^^^^ (20) 

The partial temperature and strain Tf , r = warp, weft, concentration fac- 
tors are derived from the solution of the Eshelby equivalent inclusion problem, in 
which an isolated inclusion of an ellipsoidal shape is embedded into the matrix 
and subjected to either average temperature gradient hi or strain Si found in the 
matrix. Bearing in mind certain randomness in the geometry of a single ply unit 
cell. Fig. 7(a), it is possible to derive a certain statistically equivalent ellipsoidal 
inclusion, for which the macroscopic estimates, Eqs. (18) and (20), are reasonably 
close to EE simulations for a certain range of parameters a, b, g, h. Eor the present 
material system, the three semi-axis = 1,^2 = 0.1,^3 = 0.01, characterizing 
the shape of the ellipsoidal inclusion, were found optimal regardless of the type 
of the homogenization problem. Further to this matter, searching for an optimal 
shape of the equivalent inclusion for various geometries permitted us to relate the 
values of ^2,^3 axes, given ^1 = 1, to g/a and b/a ratios, respectively, see also 
discussions in [12, 24]. Therefore, knowing at least the averages of parameters 
a, b and g allows us to define the shape of the ellipsoid with a relative ease as 

^2^^(1--), 6^^(1-4-). (21) 
la b\j a 

The subscript m in Eqs. (18), (20) identifying the mesoscopic effective prop- 
erties of a pore-free textile ply was chosen purposely as these properties are as- 
sumed to play the role of matrix in the second homogenization step, in which the 
porous phase is introduced into the new homogeneous, but orthotropic, matrix. In 
this simple case the estimates of the effective properties of the homogeneous ply 
simplify as 

X^^ - + ^2 (/tvoid ~/tm"^) Avoid' (22) 

= Lr(l-C2Af„J, (23) 
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with ^ 

Lvoid = and A^^^^ = TJ^j^ [c^I + CvoidT^oid] ' s^O.s, (24) 

Here, an optimal shape of the eUipsoidal inclusion (void) was found directly from 
2D images of the composite, see e.g. Figure 1. Notice essentially three distinct 
plies with different porosity. While the most bottom ply is free of pores (I) with 

, LJJJ^ properties, the middle (II) and top (III) ply effective properties are given 
by Eqs. (22), (23). Similar to FE simulations presented in the previous section, an 
oblate spheroid with = ^2 and ^3/^1 = 0.33, Cvoid = 0.07 for the middle ply and 
^3/^1 = 0.66, Cvoid = 0.15 for the top ply was used. 

Finally, the macroscopic properties of the homogeneous laminated plate are 
obtained applying standard averaging rules. While simple arithmetic and inverse 
rules of mixture are used for the components of effective conductivity matrix, 
the homogenized 6x6 stiffness matrix can be determined as suggested in [62, 
page 163]. Final results are summarized in Table 9 adopting the material proper- 
ties of individual phases from Tables 2-4. 



Thermal conductivities 


Elastic moduli 


[Wm-iK-i] 


[GPa] 


warp/ weft A^trans 


En G\2 £'33 Gn 


8.63 2.25 


55.7 7.8 16.3 6.5 



Table 9: Mori Tanaka estimates of the homogenized thermal and elastic moduli of the four layer 
C/C composite laminate with I-II-III-I stacking sequence as seen in Figure 1 



5. Conclusions 

The present article summarizes recent developments in the study of imperfect 
carbon-carbon textile composites initiated already in [1]. To arrive at the present 
stage of understanding the complex structural response of these material systems, 
the machinery of homogenization tools has been fully exploited revealing both 
advantages and drawbacks of individual methods. 

When introduced in the framework of hierarchical modeling, the Mori-Tanaka 
method yields reasonably accurate approximations of the effective properties yet 
at the fraction of time when compared to finite element based simulations. 

In view of real material samples with a large amount of flaws (transverse and 
delamination cracks, large intertow vacuoles) even prior to loading, the model- 
ing strategy based on the well known concept of periodic unit cell may become 
preferable, particularly if the response of the material exceeding its elastic limit 
becomes the primary interest. Presently, its potential is seen mainly in the formu- 
lation of statistically equivalent periodic unit cell that attempts to accommodate 
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the most severe imperfections of the real material. Supported by the results de- 
rived in the course of this work, the two-layer SEPUC enhanced by incorporating 
the porous phase appears to be a suitable candidate for the computational model of 
plain weave textile reinforcement based composites including, apart from the in- 
vestigated C/C composites, a large group of textile reinforced ceramics with their 
anticipated application in bio-medicine. 

To appreciate the concept of SEPUC we present, in addition, the results ob- 
tained from an independent stochastic analysis. There, a family of representative 
PUCs was generated assuming a random variation of individual geometrical pa- 
rameters of the cell including the major porosity. Each variable was assigned the 
Gaussian probability distribution with the mean values and standard deviations 
taken from Table 1. Ten such unit cells were generated using the Latin hypercube 
sampling method. Details, including the assumed correlation between individual 
variables, are available in the doctoral thesis of the first author [63]. 

All results are presented in Figure 12 confirming the predictive superiority of 
the SEPUC over other methods employed in this study. The high computational 
cost must, however, be reminded to keep the two approaches (unit cell models 
against analytical micromechanics) on the same starting line. 



12 



6 



o FEM 
* SEPUC 
MT 
A EXP 



Parameter 
(a) 



80r 



60 



§40 



20 



o FEM 
* SEPUC 
MT 
A EXP 



i>^^ C>^' iy^' C>^ 
Parameter 



(b) 



Figure 12: Comparison of the numerical and experimental results; (a) Coefficients of thermal 
conductivity, (b) elastic properties 
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